Viability of competing field theories for the driven lattice gas. 
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o ; 

D , It has recently been suggested that the driven lattice gas should be described by a novel field 

^"*^ ■ theory in the limit of infinite drive. We review the original and the new field theory, invoking several 

well-documented key features of the microscopies. Since the new field theory fails to reproduce 
these characteristics, we argue that it cannot serve as a viable description of the driven lattice gas. 
Recent results, for the critical exponents associated with this theory, are re-analyzed and shown to 
be incorrect. 
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The critical behavior of the driven lattice gas (DLG) |l]] has been the subject of some debate, ever since the first 
Monte Carlo simulations Q and field theoretic predictions were found to give differing values for the order 
parameter exponent p. This discrepancy has led to developments in different directions: some researchers have 
modified the simulation data analysis, invoking anisotropic finite size scaling 0, while others H|| have suggested 
that the original field theory might be deficient in the limit of infinite drive, proposing |||| and analyzing JTo| an 
alternate coarse-grained theory instead. In this communication, we review both the original |}|4j and the alternate 
|§|-|lof field theory, in the light of Monte Carlo simulation data. We first document that the alternate theory is not a 
I ' coarse-grained description of the driven lattice gas, since it fails to exhibit several well-established properties of the 
""O • microscopic model. In a second step, we re-analyze the proposed theory, assuming that it might describe some other, 
yet to be determined, microscopies. We show that the renormalization group analysis of Ref. |l0| ] is seriously flawed, 
resulting in incorrect exponents and a proliferation of uncontrolled infrared singularities. 

We begin with a brief summary of the background. Microscopically, the DLG is a simple ferromagnetic Ising lattice 
gas, half-filled and coupled to a heat bath at temperature T, in which particles jump to empty nearest-neighbor 
sites subject to the usual Ising energetics and a uniform driving force E acting along a particular lattice direction. 
Thus, the effect of E is identical to adding a locally linear potential. Clearly, E — corresponds to the equilibrium 
Ising model. On the other hand, even E = oo can be realized if Metropolis rates are used: Simply accept/forbid 
^v^j ' all forward/backwards jumps. Since large values of E accentuate the nonequilibrium features of this system, most 
£vq , simulations have been performed at E > 50, in units of the Ising coupling constant. 

t-H ' The driven lattice gas and many of its variants have attracted considerable attention since they evolve into simple 
0^ \ nonequilibrium steady states displaying a wealth of counterintuitive characteristics jllj . Two of its most remarkable 
features are (i) the discontinuity singularity of the structure factor S(k) p],|l2|], which is intimately connected to an 
r~ d decay (in d dimensions) of the two-point correlations Jl^Jbfl , and (ii) the emergence of nontrivial three-point 
correlations in the disordered phase, corresponding to the violation of the Ising symmetry by E (which drives 
particles and holes in opposite directions) . Such dramatically "non-Ising" characteristics are easily observed in Monte 
Carlo simulations, at intermediate and large driving fields. They are also confirmed in a high-temperature series 
expansion, derived directly from the microscopic dynamics |l3|]l6f |. 

These observations from Monte Carlo simulations play a crucial role in identifying the correct field theory. A basic 
O ■ tenet in the study of critical phenomena is that a microscopic model and its coarse-grained field theory should possess 
the same symmetries, if they are to belong into the same universality class. For the driven lattice gas, the data on the 
structure factor indicate that the theory is highly anisotropic. Moreover, the detailed behavior of the discontinuity 
singularity, upon approaching the origin in wave vector space from different directions, informs us precisely how the 
familiar Ornstein-Zernike form is modified. Generically, we find pd[ | that 

R = ^IkLHoSQcj-.*!! =°) > 1 
limfc^o S(k± = 0,k\\) 

above criticality, and R — > oo upon approaching T c . The subscripts distinguish the parallel (||) and transverse 
(_L) subspaces, with respect to the drive direction. Just as significantly, the non-vanishing three-point functions 
demonstrate that the usual "up-down" symmetry of the Ising model is broken. 

These key features of the microscopies must be reflected in any viable continuum theory for the driven lattice gas. 
We first consider the original field theory It is based on a Langevin equation, in continuous space and time, 
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which describes the stochastic evolution of the local particle density p(x, t). In terms of 4> = 2p — 1, the equation 
reads: 



^0 = A{(r ± -Vi)Vi^ + T||Vf0 + 5V||0 2 + |vi ( /) 3 }-Ve . (2) 
The Langevin noise term reflects the fast degrees of freedom: 

(£(x,t)> =0 

(V£(x,t)V'£(x',t')) =-2(n±V 2 ± +n ll V^6{x-x')6(t-t') . 

We emphasize that (i) all coefficients are strictly positive except possibly t± and/or T|| which control criticality (see 
below) and (m) independent from one another (i.e., not related by symmetry). The parameter A sets the time scale. 

This theory contains two closely linked key ingredients: First, there is a driving term, £V||</> 2 , where £ denotes 
the coarse-grained drive (a naive continuum limit gives £ oc tanh(£'/X')). This term is required to break the Ising 
"up-down" (<ft — > —<ft) symmetry. Second, the theory is highly anisotropic, with two different diffusion coefficients T|| 
and tj_. In particular, it predicts an equal-time structure factor, 

+ fc 2 

^ (k) = r ± kl +r||fc 2 + 0(fc4) ^ 

in the disordered phase. This S generates a discontinuity singularity R = (n±r^ ) / (n^Tj_). To ensure that the observed 
behavior is faithfully reproduced, we demand nj_r^ > ni\T± in the disordered phase. Moreover, criticality must be 
marked by tj_ = at positive t\\ if the divergence of 5* is to be captured correctly. To summarize, the two key features 
of the original Langevin equation are unambiguously supported by the Monte Carlo data for the microscopic model. 

We comment briefly on the issue of finite versus infinite fields. In all Monte Carlo simulations, the current is 
observed to saturate as E increases. This saturation is reflected by £ oc tanh(£'/T) — > 1 in the original field theory. 
Therefore, this theory holds equally well for any nonzero value of the microscopic drive. Furthermore, simulations 
using Metropolis rates with E = 50, 100 and oo have been performed. The results are (statistically) indistinguishable! 
Such sensible behavior is entirely consistent with this theory. 

The discrepancies arise when critical exponents are measured, specifically the order parameter exponent /3, and 
compared to field theoretic predictions. The original field theory, due to the vanishing of t_l at positive tii , naturally 
leads to anisotropic scaling of wave vectors: k\\ ~ fc^ +A in the critical region, with a nontrivial anisotropy exponent A. 
Three important consequences are that, first, the upper critical dimension d c is shifted to 5, and second, the theory 
predicts A = 1 + (5 — d)/3 and (3 = 1/2 exactly, i.e., to all orders in perturbation theory. The values obtained by 
simulations differ, depending on the method used to analyze the data. If a careful anisotropic finite size analysis M 
is used, based on system sizes consistent with the field-theoretic scaling, i.e., i||/i 1 L +A = const, the field-theoretic 
exponents result in good data collapse for a number of different observables ^,0. However, data for isotropic systems, 
Lu/L± = const, appear to indicate an order parameter exponent around 0.23 |2|. Since most of the data were taken 
at very large fields, some authors |^|| have suggested that the origin of the discrepancies does not reside in the data 
analysis. Instead, they argue that the standard field theory does not capture the E — > oo limit correctly and propose 
an alternate theory. It is based on the Langevin equation: 

^0 = A{(r ± -Vi)Vi0-VfVi0+|vi ( /» 3 }-Ve . (4) 

With minor renamings of parameters |17| , this is Eq. (1) of Ref. fl0|| . The vanishing of t± marks the critical point. 
The noise satisfies (Eq. (2) of Ref. 0): 

<£(x,<)) = o 

<V£(x, t)V'£(x', t')) = -2A (Vi + iv 2 ) <5(x - x ')S(t - t') . 



Two key terms appearing in the original field theory are absent in this one, namely, 

• the driving term V||</> 2 , and 

• a diffusion term T|| V 2 </> for the parallel direction. 
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Since the driving term is absent, the alternate field theory obeys the Ising "up-down" (<f> — > —cf>) symmetry. Thus, 
three-point functions are identically zero in this theory, for all T >T C . This prediction is in serious disagreement with 
existing Monte Carlo data! While one may argue that a field theory need not reproduce all of the microscopic detail 
of the underlying lattice model, one should be very cautious before endowing it with a higher symmetry: This is only 
justified if a high-symmetry fixed point exists and can be shown, via an explicit renormalization group calculation, to 
be stable against perturbations by symmetry-breaking operators. Neither is the case here. 

The absence of the parallel diffusion term also has serious consequences. Eq. (|]) generates a steady-state structure 
factor: 



s(k) = 



ki 



1 7,2 

2 K || 



k\ (k 2 + rj.) 



(5) 



which ought to be a good approximation at high temperatures. Yet, for ku ^ it predicts a divergence along the 
whole k± = line, at any T > T c . This stands in glaring contrast to the Monte Carlo results for the disordered phase, 
where all measured structure factors are found to be finite everywhere in fc-space. 

Since Eq. (|]) fails to reproduce the most basic properties of the microscopic model, we conclude that it is not a 
viable field theory for the driven lattice gas. It may, however, describe some as yet unknown microscopies. Therefore, 
we now proceed to analyze the field theory, defined by Eq. (^), in its own right. 

Following Ref . M , we recast Eq. ([|) as a dynamic functional JT^] : 



d d xdt 



d t ci> + \{yW 2 



_L 



Xcj) V 



ft 



(6) 



We first note that Eq. (||) describes a theory with a four-point coupling 0V^0 3 and anisotropic free propagators 
as given in Ref. Therefore, the combinatorics of this theory is identical to that of Model B, which reduces to 

(/> 4 -theory in the static limit. For such theories, it is well known fli| that the one-loop result for the exponent v 
(denoted u± in Ref. ]l0| ) is determined by combinatorics alone, i.e., the explicit expressions for the Feynman integrals 
are not required. This is most easily seen by calculating in the critical theory, where tj_ = 0, with insertions of 
A0Vj_<^. We denote one-point irreducible vertex functions with n (n) external <f> (</>) legs and m insertions by r~JJ . 

At one-loop order, there are two primitively divergent vertex functions, namely and . Both of these consist 
of a zero-loop term and a one-loop contribution. Each one-loop contribution consists of a combinatoric factor, the 
appropriate powers of the coupling constant and the external momentum, and a loop integral. The key simplification 
here is that the loop integrals for and are identical, independent of the detailed forms of the free propagators. 
Thus, the two one-loop contributions differ only by a simple factor which is purely combinatoric in origin. As a result, 
one obtains to first order in e = d r — d, for all of these theories: 



(7) 



Since the authors of Ref. |To| have chosen to calculate at finite t± , let us illustrate how this result emerges in their 
case. No insertions are needed here so the upper index of can be dropped. Keeping track of coupling constants 
and signs, and taking care of the T c shift, we can write the two bare vertex functions fn and Ti3 in the form 



Tn = iuj + Xk\k 2 + Xk z ,r± 



1 + \gh 



Ti 3 = Xgk A 



(8) 



Here, the factors 1/2 and —3/2 arise from combinatorics while the integrals Ii and I 2 are easily computed in dimen- 
sional regularization, resulting in 



h 



(47T) 



[l + 0{e)\ and J 2 



(47T) 



[l + 0(e)] 



(9) 



We notice immediately that the simple e-poles of I\ and I2 are identical. Thus, their numerical prefactor can be 
absorbed into the definition of the coupling constant, leaving us with one-loop corrections to v that are purely 
combinatoric in origin. Completing the calculation at finite t±, this provides the key to Eq. (J?]). Only at two-loop 
order do the detailed forms of the free propagators come into play. Then, of course, exponents are also no longer 
determined by combinatorics alone. 
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In Ref. (T(]], the exponent v is quoted as (1 + e/4)/2, indicating the presence of a computational error. More 
seriously, however, there are deeper flaws in this theory. Recall that the steady-state structure factor, Eq. (||), 
diverges along the whole k± = line, even for tj_ > 0. As a result, the theory is plagued by infrared singularities, 
which are entirely unrelated to criticality, and unrenormalizable divergences. We note, for completeness, that one can, 
of course, regularize such singularities by re-introducing the diffusion term T|| Vm<?!> into Eq. Then, however, one 
should also reconsider the two fourth-order derivative terms. At two-loop order, these will acquire different primitive 
divergences, so that an additional coupling constant p 2 is required, appearing in Eq. (||) as p 2 (V^) 2 + V|Vj_0. 

To summarize, we have shown that the field theory proposed by Garrido, de los Santos and Muhoz |||| fails 
to reproduce the key features of the driven lattice gas. Predicting infinite structure factors and zero three point 
correlations (for all temperatures above criticality), it cannot be a viable continuum model for the latter. Accepting 
it as a representation of some other, as yet undetermined, microscopic model, we carry out a standard analysis. 
First, we find that the one-loop calculation of Ref. JhJ is incorrect. Second, beyond one-loop order, uncontrolled 
infrared singularities proliferate, rendering the field theory unrenormalizable. In contrast, the original field theory 
HQ is consistent with the fundamental symmetries of the driven lattice gas, for any value of the drive. Its predictions 
for 2- and 3-point functions in the disordered phase are in good agreement with simulation results. Based on the 
phenomenology of S(k) near criticality, it plumbs the consequences of a highly anisotropic scaling limit, fc|| ~ k^ A — > 
0. To test its predictions against Monte Carlo simulations, this limit should be respected in the choice of system 
sizes, i.e., L« ~ L^ +A . If, instead, simulations and finite-size analysis are performed with disregard for such strong 
anisotropies, complications from extraneous scaling variables Q or inconsistencies |2C| can be expected. In a more 
exotic scenario, such simulations may indicate a new type of low temperature phase, quite distinct from the ordinary 
Ising-like one. 

This work is supported in part by the National Science Foundation through the Division of Materials Research. 



[1] S. Katz, J. L. Lebowitz, and H. Spohn, Phys. Rev. B 28, 1655 (1983); J. Stat. Phys. 34, 497 (1984). 

[2] J. Marro, J.L. Varies, and J.M. Gonzalez-Miranda, Phys. Rev. B 35, 3372 (1987); J.L. Varies and J. Marro, J. Stat. Phys. 

49, 89 (1987); J.-S. Wang, K. Binder, and J.L. Lebowitz, J. Stat. Phys. 56, 783 (1989). 
[3] H.K. Janssen and B. Schmittmann, Z. Phys. B64, 503 (1986). 

[4] K-t. Leung and J.L. Cardy, J. Stat. Phys. 44, 567 and 45, 1087 (Erratum) (1986). 
[5] K-t. Leung, Phys. Rev. Lett. 66, 453 (1991) and Int. J. Mod. Phys. C3, 367 (1992). 
[6] J.-S. Wang, J. Stat. Phys. 82, 1409 (1996). 

[7] K. Binder, in: Finite Size Scaling and Numerical Simulations of Statistical Systems, ed. V. Privman (World Scientific, 
Singapore, 1990). 

[8] P.L. Garrido, F. de los Santos, and M.A. Munoz, Phys. Rev. E 57, 752 (1998) 
[9] F. de los Santos and P.L. Garrido, J. Stat. Phys. 96, 303 (1999). 
[10] F. de los Santos and M.A. Munoz, to appear in Phys. Rev. E 61 (February 2000), |cond-mat/9910295| (1999). 
[11] For a more extensive review, see, e.g., B. Schmittmann and R. K. P. Zia, Phase Transitions and Critical Phenomena, Vol. 

17, edited by C. Domb and J. L. Lebowitz (Academic, London, 1995). 
[12] R.K.P. Zia, K. Hwang, K-t. Leung, and B. Schmittmann, in Computer Simulations Studies in Condensed Matter Physics 

V, eds. D.P. Landau, K.K. Mon, and H.-B. Schiittler (Springer, Berlin 1993); see also Fig. 3.2 in Ref. |n). 
[13] M.Q. Zhang, J.-S. Wang, J.L. Lebowitz, and J.L. Valles, J. Stat. Phys. 52, 1461 (1988). 
[14] P.L. Garrido, J.L. Lebowitz, C. Maes, and H. Spohn, Phys. Rev. A 42, 1954 (1990). 

[15] K. Hwang, B. Schmittmann, and R. K. P. Zia, Phys. Rev. Lett. 67, 326 (1991) and Phys. Rev. E 48, 800 (1993). 
[16] B. Schmittmann and R. K. P. Zia, J. Stat. Phys. 91, 525 (1998). 

[17] In Ref. [px|, A is denoted by e /2 and r± is just r. Our renamings make comparison of Eqs. (|) and M) easier. 

[18] H.K. Janssen, Z. Phys. B: Cond. Mat. 23, 377 (1976); C. De Dominicis, J. Phys. (France) Colloq. 37, C247 (1976); 

H.K. Janssen, in From Phase Transitions to Chaos, G. Gyorgyi, I. Kondor, L. Sasvari, and T. Tel, eds. (World Scientific, 

Singapore, 1992). 

[19] See, e.g., D.J. Amit, Field Theory, the Renormalization Group and Critical Phenomena, 2nd revised edition, (World 
Scientific, Singapore, 1984); J. Zinn- Justin, Quantum Field Theory and Critical Phenomena, (Oxford University Press, 
Oxford, 1989). 

[20] K.-t. Leung and J.S. Wang, Int. J. Mod. Phys. C 10, 853 (1999). 



4 



